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. . . Abstract 

, We calculate the energy-momentum tensor due to electromagnetic vacuum fluc- 

■ tuations between two parallel hypcrplanes in more than four dimensions, considering 
both metallic and MIT boundary conditions. Using the axial gauge, the problem can 
be mapped upon the corresponding problem with a massless, scalar field satisfying 
respectively Dirichlet or Neumann boundary conditions. The pressure between the 
plates is constant while the energy density is found to diverge at the boundaries when 
there are extra dimensions. This can be related to the fact that Maxwell theory is 
then no longer conformally invariant. A similar behavior is known for the scalar field 
where a constant energy density consistent with the pressure can be obtained by im- 
proving the energy-momentum tensor with the Huggins term. This is not possible for 
the Maxwell field. However, the change in the energy-momentum tensor with distance 

■ between boundaries is finite in all cases. 

O 

' 1 Introduction 

When a classical field is quantized, the modes that can be excited are solutions of the 
classical wave equation and are labeled by different quantum numbers. These modes will 
depend on the imposed boundary conditions and will therefore be influenced by the pres- 
ence of confining boundaries. Each such mode has a zero-point energy which contributes 
to the total vacuum energy of the field. As a result, these vacuum fluctuations give the 
ground state of the system an energy which depends on the presence of nearby boundaries. 
For two parallel and perfectly conducting plates placed in vacuum CasimirpQ showed that 
for the electromagnetic field this energy corresponds to an attractive pressure 

P = -— (1) 
240L 4 { J 

between the plates, separated by a distance L. This macroscopic quantum effect was 
for a long time in doubt, even after the first experimental verifications by Sparnaayp]. 
Today this Casimir force is measured to high precision [3] and even effects of non-zero 
temperatures are being investigated [1]. 
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Together with this experimental progress, modern regularization methods to remove the 
unphysical divergences endemic in these calculations now make them much simpler than 
previously [5J. One can then investigate more detailed properties of the effect like how 
the energy or stresses are distributed between the plates. One must then calculate the 
vacuum expectation value of the full energy-momentum tensor as was done by Liitken and 
Ravndal[6]. This confirmed a previous calculation by de Witt [7] of the fluctuations of the 
electric and magnetic fields near a metallic boundary where they diverge but in such a way 
that the energy density remains finite. Since Maxwell theory is conformally invariant in 
D = 4 spacetime dimensions, one can directly relate this energy density to the attractive 
pressure [8]. 

Recently the Casimir energy has been invoked to explain the dark energy which seems to 
drive the present acceleration of the Universe [9]. It appears in particular in cosmological 
models with extra dimensions [10]. In these models the confinement of the fluctuating 
fields is provided by compactification of the extra dimensions. The energy appears as 
a cosmological constant in our four-dimensional Universe and is given numerically by a 
formula of the same form as ([1]) with L given by the size of the compactified dimensions. 

The calculation of Casimir energies in higher-dimensional spacetimes was first done by 
Ambj0rn and Wolfram|llj. They calculated the global energies equivalent to the Casimir 
force and thus obtained no knowledge of how the energy is distributed between the hyper- 
planes. Since the force due to electromagnetic fluctuations is expected to be proportional 
to the force due to fluctuations of a massless scalar field, only the effects of this kinemat- 
ically simpler field was investigated. Each mode has a momentum transverse to the 
normal of the plates plus a component k z = nir/L with n = 1, 2, . . . in the direction of the 
normal. If d is the number of spatial dimensions, the vacuum energy between the plates 
per (d — l)-dimensional hyperarea will follow from the divergent integral 



We can now do the transverse integration by dimensional regularization. The remaining 
sum over n is then done by analytical continuation of the Riemann zeta function Cr( z ) to 
give 



for the zeta function. Then we can write E = SqL when we introduce the energy density 






(3) 



It can be simplified using the reflection formula 




(4) 



F(D/2)Cr(D) 



(5) 
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Defining now the pressure between the plates by P = —dE/dL, it is then simply 

P = (D - 1)£ (6) 

where D = 1 + d is the spacetime dimension. Taking D = 4 and multiplying the result by 
two for the two polarization degrees of the photon, we recover the original result ([1]). In 
D-dimensional spacetime we must for the same reason multiply the result ([6]) by D — 2 to 
obtain the electromagnetic Casimir force. 

The spatial distribution of the vacuum energy can be obtained from the energy-momentum 
tensor. For the scalar field it is 

Tp.v = d^4>dv4> - -q^C (7) 

where the massless Lagrangian is C = (l/2)(d\(ft) 2 choosing the metric to be rjuv = 
diag(l, -1, • • • ,-1). Since its trace = (1 - D/2)(d X (p) 2 is zero only in D = 2 di- 
mensions, it is in general not conformally invariant in higher dimensions. Calculating now 
the Casimir energy density (Too) i n f° r example D = 4 spacetime dimensions, one then 
obtains a result which diverges at the plates after regularization[T2]. When integrated, it 
will thus not reproduce the total Casimir energy corresponding to the force ([6]). 

For the scalar field this apparent problem can be solved. One can improve [13] the above 
energy-momentum tensor in any spacetime dimension D > 2 by adding the Huggins 
termpH] 

AT^ = -I:^:^ - v , u d 2 )4> 2 (8) 

The energy-momentum tensor is then traceless and conformal invariance has been restored. 
Since this new term is a divergence, it will not contribute to the force, but change the 
distribution of energy around the plates. Vacuum expectation values of all components 
of the energy-momentum tensor are now constant between the plates and zero outside as 
already noticed by de Witt[7j, Milton [15] and others [12]. 

However, the situation for the electromagnetic field is somewhat different. It has the 
energy-momentum tensor 

= F m F a v + -fl^F^F* (9) 

where F^ = d^A u — d u A^ is the Faraday tensor. From the trace T M M = (— 1 + D/A)F 2 j3 
we see that it is conformally invariant only in D = 4 spacetime dimensions. In this case 
the Casimir energy is also constant between the plates[H] as for the scalar field. But for 
dimensions D > 4 it is no longer clear how the energy is distributed since there is no way 
to construct a gauge-invariant analogue of the Huggins term in Maxwell theory. 

In the following we will investigate this problem in more detail. In D = 4 dimensions one 
can choose the transverse gauge and expand the classical field in electromagnetic multi- 
poles. This is cumbersome in higher dimensions since the field then has more magnetic 
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than electric components. Prom the geometry of the problem it is more natural to choose 
the axial gauge n^A^ = A z = where the D-vector = (0, 0, . . . , 0, 1) is normal to the 
plates and is called the z-direction. The Faraday tensor then has a correlator which can 
be directly obtained from scalar field theory in the same geometry. In the next chapter 
we therefore derive a general expression for the scalar field correlators, both in the case 
of Neumann and Dirichlet boundary conditions. These can then be used to calculate 
the expectation value of the scalar energy-momentum tensor (|7J) and the relevance of the 
Huggins term is discussed. 

For the electromagnetic field considered in Chapter 3, we need to know the boundary 
conditions. These can be of the metallic type used for the standard Casimir force in 
D = 4 or the QCD version used in the MIT bag model for confinement of quarks. In 
the axial gauge we find that these two possibilties correspond to Dirichlet and Neumann 
conditions for the corresponding scalar field. The fluctuations of the different components 
of the Faraday tensor are then calculated with particular attention to the energy density 
and the pressure between the plates. While the pressure is found to be constant and in 
agreement with the global Casimir force, the energy density diverges at the plates. 

In the last chapter this problem, which no longer can be cured with a Huggins term, is dis- 
cussed and compared with similar divergences in other systems. With physical boundaries 
that only confines fluctuations with frequencies below a certain cut-off, all field fluctuations 
should reach a finite value when the boundaries are approached. 



2 Scalar fields 



It will be very convenient to denote by a bar any vector or tensor component orthogonal 
to the unit normal of the plates, which is taken to be in the z-direction. Thus for a 
full D-vector we write A = (A») = (A* 1 , A z ). Note that the (D - l)-vector A = (A^) also 
includes the time component A . The metric can thus be written as if"* = fj 



where fj is the projection of rj onto the barred subspace. The field operator for a massless 
scalar field satisfying the Dirichlet boundary conditions (f>(x, z = 0, L) = 0inD = d+ l 
spacetime dimensions will then be 



x 



\fl E / ^^\f^ M ky ) e ~ fk '" + ^) Mnnz/L) (10) 



Here k = (uj n , ky) with the frequency u n = y ky + where m n = k z = mr/L. Had we 
instead chosen the Neumann boundary condition d z (p(x, z = 0, L) =0, we just have to 
make the replacement sm(nirz/L) — > cos(nirz/L) in the sum. Then there should also be 
a n = mode to be included in the sum with normalization constant But with 

the regularization we will use in the following, it will not contribute and is therefore not 
further considered. 
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2.1 Feynman correlator 



In the above field operator for the scalar Dirichlet modes we have used a normalization 
which corresponds to [a n Qs.T)> oL(k^)] = 5 nn i{2-n) d ~ 1 5i}^T — ky) for the annihilation and 
creation operators. From this we find the Feynman propagator 

G D (x,x') = {VL D \T4>{x)4>{x')\n D ) 

d d k 2 ^ sin(nnz/L) sm(mrz'/L) (u) 
(2ir) d L^ k 2 -ml + ie 1 ' 

n=l " 

Here we integrate over all components of the tf-dimensional Lorentz vector k. Assuming 
e = x — x' to be spacelike, we may choose a coordinate system where it has no components 
in the time direction. We can then rotate k° to the imaginary axis and find 

, 2 f d d k ^ sm(nirz/L)sm(mrz'/L) ^. e - 

Gd{x > x ] = ij (2^^ wT^ n 6 (12) 

(where now A; is a Euclidean vector. In order to evaluate the sum, we consider the function 

2 sm(nirz/L) sm(nirz'/L) 

n=l 

which solves the differential equation 

~ + lAg(z 1 zf) = S(z-z') (14) 

on [0, L] with Dirichlet boundary conditions. As can be verified by insertion, the solution 
is 

sinh kz < sinh k(L — z>) 

9{z, z ) = . - (15) 

k sinn kL 

where z < = min(z, z') and z> = max(z, z'). Expanding the hyperbolic functions, one then 
finds 



. Z, Sill 71/1 ^ / Lj I 51111 U ll 6 1^ 

^) = z^ w^i (13) 



1 oo 

gfa z 'j = _L L-k\z-z'\ _ e -fc(«+a') _ e -fe(2L- z - 2 ') + e -fc(2L-| Z - Z '|)\ ^ g -2 



i=o 



^ / 27T fc2 + fc 2 ^ 



j=-oo 



The last equality is verified by evaluating the k z integral by contour integration. Inserting 
now this partial result into (|12p and using rotational invariance, we find 

G D (x,x') = E / (OT^F (e^- z '- 2 ^+^ - e i ^^+-'- 2 ^)+^-]) (17) 



j=-oo 
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where now the (d + l)-dimensional vector k = (k, k z ). Each of the integrals are given by 
the generalized Coulomb potential 

_ f d n k 1 , k . (x „ x >) _ r(n/2 - 1) 



T4(x x ) - J (27r)n ^ e J ■ x x - 4W2|x _ x/|re _ 2 (18) 

in n = d+1 spatial dimensions. Had we instead considered Neumann boundary conditions, 
the sine function in (jlOf) would have been replaced by the corresponding cosine function. 
The only change would then have been that the last term in (|17p came in with opposite 
sign. Introducing the D-vectors Zj = (e,z — z' — 2jL) and zj = (e, z + z' — 2jL) of lengths 
Rj = {z 2 -) 1 ! 2 and Rj = (zj) 1 / 2 , we can now write the result for both correlators as 

oo 

G(x,x') N/D = [Vd+i(Rj) ±V dn (R j )] (19) 

j=-oo 

where the upper sign is for Neumann and the lower for Dirichlet boundary conditions. 
For a massive field we would have found a similar result, but with the Coulomb potential 
replaced by the corresponding generalized Yukawa potential. 

In fact, almost every student of introductory electrostatics could have written down this 
result immediately by realizing that the problem is equivalent to calculating the potential 
of a point charge between parallel plates in D = d+1 spatial dimensions, using the method 
of images to enforce the boundary conditions. 



2.2 Energy-momentum tensor 

The term V^+i(i?o) hi (fTU|) is equal to the free correlator Gq(x,x') = {0\Tcj)(x)4>(x')\0) , 
where 1 0) is the bulk vacuum. It diverges in the limit x' — > x. But defining now the 
physical vacuum expectation value of the energy-momentum tensor by the point-split 
limit 

(T^{x)) = lim[(fi| Tfja/ (x , x - (0\T^(x',x)\0)] (20) 

x' — *x 

its contribution is removed. The finite expectation values will then follow from the reg- 
ularized correlator Gd(x,x') — Gq(x,x') which contains the effects of the plates. We will 
continue to denote it by Gd(x, x') in the following and it is given by (|19p when we in the 
first part leave out the j = term. Since we have assumed that x — x' is non-zero and 
spacelike, the time-ordering symbol in the correlator can be ignored and it satisfies the 
Klein-Gordon equation (B 2 — d 2 + m 2 )G(x, x') = for both boundary conditions. 

For the Dirichlet vacuum expectation value of the scalar energy momentum tensor ([7]) , we 
first need the part 

{d%<j){x)dl,<i){x')) D = -aga%G D (x,x?) = -^ d B 2 g d {x,x') 

= \rf v {m 2 -d 2 z )G D {x,x') (21) 
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using Lorentz invariance. Similarly, it follows that 

( d z 0(x)d z ,<P(x')) D = -8 2 G N {x, x') (22) 



since the two parts in the correlator (|19|) has opposite symmetry under the exchange 
z — > z'. For vacuum expectation value of the point-split Lagrangian 

C(x,x') = i[^^(x)fl^0(a/) - d z <t>{x)d z ,<t>{x') - m 2 <P(x)<P(x')} (23) 
we thus find 

(C{x,x')) D = ^d 2 z [G N (x,x')-G D (x,x')\ (24) 

The point-split expressions for the canonical energy-momentum tensor ([7]) are thus found 
to be 

(7» D = ^(^Gd-^Gd + ^Gn-Gd)]) (25) 

(T ZZ ) D = -±%(G N + G D ) (26) 

Corresponding results for the Neumann expectation values are obtained by the exchange 
D <-> N. The physical limit x — > x' can now be taken where a resulting z-dependence 
can only come from the last sum in the correlators (|19p . But for the pressure P = (T zz ) 
we see that this will cancel out in the sum Gd + Gfq so that the pressure is constant 
between the plates. This is physical reasonable and is also the case for the fluctuations of 
a massive field. It follows directly from the conservation of the energy-momentum tensor. 
The expectation values of the other other components of the energy-momentum tensor in 
(|25p will in general be dependent on the position z between the plates. 



Let us now calculate the pressure in the massless limit. We will then need the double 
derivative 8 2 {Gd + Gn) which follows directly from (|19l) in the limit X ^ X clS 

dl(G N + Gd) = 2 hm d 2 V d+1 ( Rj ) = 2d(d - 1) ^'Jj £ 

j=—oo j=—oo 

2{D - l)T(D/2) 

~ WI*IP UD) (2?) 

where the ' denotes that j = is excluded from sum. Using this in (|26|) we reproduce 
exactly the standard pressure © obtained from the total energy. It is seen to be the same 
for both boundary conditions. 

The energy density 6 = { Too) between the plates follows from (|25|) . When the mass m = 0, 
we then need to calculate in addition the quantity 

oo 



d 2 z (G N -G D ) = 2Km ^ d^Vd+xiRj 



j-- 



(D-l)T(D/2) 
(47r) D/2 L D fo{z/L) (28) 



7 



when we introduce the function 

oo 1 

= E FTW <29) 

3=-oo 

Notice that the term j = is now to be included. The sum can be expressed by the 
Hurwitz zeta function 

( H (s,a)=J2 - ( T u ( 3 °) 

which allows us to write 



M^) = Ettttttw + E 



- (j + z/L) D ^ (j - 
Ch(D,z/L) + Ch(D,1-z/L) (31) 



When Z? = 4 the same position-dependent term was derived on this form by Kimball[15j. 
But when the spacetime dimension D is an even number, we can express the result in 
terms of the digamma function ip(x) using the relation 

(-l) k ( d \ k ~ 1 



We then have 

f D (z/L) = ]_ [j- \ [^(x)-^(l-x)J ( J D = even) (33) 

where x = z/L. This simplifies even more since ip(x) — — x) = — 7rcot(7rx), which 
allows us to write 

/D(Z/L) = T(D) V COt ^ (D = even) (34) 

with = irz/L. For the ordinary Casimir effect in D = 4 spacetime dimensions, this 
function also appeared in the calculation of the electromagnetic field using another regu- 
larization and choice of gauge[6]. This follows from writing f^(z/L) = (7r 4 /3)F(6) which 
gives 

F(9) = (-^ ~ -ir) (35) 
V sin t> sin J 

This function will then characterize all position-dependent expectation values when D = 4. 

Collecting the above results, we now have the scalar vacuum energy density in arbitrary 
spacetime dimensions 



£d/n 



rp/2) 

( 47r )D/2 L D 



Cr(D) ± (D/2 - l)f D (z/L) (36) 
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where the lower sign is for Neumann boundary conditions. While the first term corresponds 
to a constant density, the last term gives a position-dependent contribution which in 
general diverges at the position of the plates, i.e. where z = and z = L. Only in the 
special case D = 2 when the scalar field has conformal invariance, will it be absent. When 
integrated over the volume between the plates, the first term alone is seen to give the 
total Casimir energy ([5]). The last term gives a divergent contribution to the same energy 
and should be absent. No such term was found in the calculation of the pressure. It 
is consistent with just the first part of the energy density which alone gives the correct 
Casimir force. 



2.3 Huggins term 

A free, massless scalar field can couple to gravity in a conformally invariant way. The re- 
sulting energy-momentum tensor will thus be tracelessfT] |13j. It differs from the canonical 
expression (|7|) by the extra Huggins term (JBJ). Using the equation of motion c? 2 = 0, one 
then finds that the improved energy- momentum tensor indeed is traceless. 

When we now want to evaluate the Huggins term for the vacuum between the two plates 
using the above point-split regularization, we interpret 

where the primed derivatives are with respect to x'. This gives the ground state expecta- 
tion value 

( AT„ U ) D = \% v ^1 d\(G N - G D ) (37) 

which is seen to be proportional to (£). The Huggins correction has no components in 
the ^-direction, leaving the pressure unaltered. The z-dependent terms cancel against the 
same terms in the canonical part (|25|) so that the resulting energy density will be constant. 
In fact, when m = 0we have for both sets of boundary conditions that 

( T„ v + ATp) = £ (V + (D- l)n M n„) (38) 

when expressed in terms of the energy density (J5J). The pressure is simply D — 1 times 
this constant energy density. This is a direct consequence of the energy-momentum tensor 
now being traceless. 



3 Maxwell fields 



The electromagnetic Lagrangian density C = —F^ u /4 is gauge invariant. In a general 



spacetime it can be written as 



£ = \e} - \b% (39) 
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where E± = — (Ai + diA°) are the components of the electric field vector while the magnetic 
field is given by the antisymmetric tensor Bij = diAj — djAi. 

We want to calculate the vacuum expectation value of the corresponding energy-moment- 
um tensor ([U]) between the plates. As already explained in the introduction, it is then 
most natural to work in the axial gauge n^A^ = 0. Since the plate normal vector n M only 
has a component along the z-axis, this requires the component A z = 0. The component 
Aq is no longer a canonical variable, but depends on the others via the Maxwell equation 
diF i0 = 0, which gives 

Aq = -A^diAi (40) 

where the operator A = df. There are thus D — 2 independent degrees of freedom in a 
D-dimensional spacetime described by the spatial field components A{ where % ^ z. We 
can then express the full Lagrangian in terms of these fields. After partial integrating and 
neglecting surface terms, we find it to be 

L = ^jd d x [A (<% - ftA" 1 ^-) Aj - Ai (did, - <%A) Aj\ (41) 

As usual, the first or electric part acts like a kinetic energy while the magnetic part acts 
like a potential energy. 



3.1 Boundary conditions and the correlator 

In order to quantize this theory, we must impose boundary conditions for the electromag- 
netic field components on the confining plates. For the original Casimir effect in d = 3 
dimensions one had metallic plates in mind where one could impose the standard con- 
straints n x E = and n • B = on the elecric and magnetic field vectors. For the more 
abstract case we have in mind here, we could just as well consider the MIT boundary con- 
dition n^F^y = proposed at MIT for the quark bag model ensuring color confinement [16] . 
In terms of components, this is equivalent to the two conditions n x B = and n • E = 0. 
They are therefore just the electromagnetic duals of the metallic boundary conditions. 

For our problem under consideration the MIT boundary condition can most directly be 
imposed. With the normal vector along the z-axis, it is equivalent to Fp z = in our 
previous index notation. Now in the axial gauge this is simply equivalent to the Neu- 
mann boundary condition d z A{(x, z = 0, L) =0, and defining A^ 1 in (|40p with Neumann 
boundary conditions. 

With this index notation, the metallic boundary conditions can also easily be generalized 
to higher dimensions by noticing that in d = 3 dimensions they correspond to Fno = 0. 
In the axial gauge this is achieved in all dimensions by requiring Ai{x, z = 0, L) = 0, i.e. 
Dirichlet conditions, and defining A -1 in (|40p with Dirichlet boundary conditions. In this 
way we can take directly over many of the previous results for the scalar field. 
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The field components Ap obey the classical wave equation d 2 Ap — dp{d v Ay) = 0. Solutions 
will be of the form Ap(x) = an{x)b(z) where the function b{z) oc cos(nirz/L) when we 
impose MIT, i.e. Neumann, boundary conditions and b(z) oc sm(mrz/L) for metallic or 
Dirichlet boundary conditions. With these boundary conditions the remaining functions 
ap(x) satisfy (d 2 + rrifyap — dp(d u ap) = which are just the equations of motion for a 
massive vector field with mass m n = nn/L. With this observation, we then immediately 
have the correlator 

D fiP (x,xf) = (n D \TAp(x)A p (x')\n D ) 

f d d k 2 ^ sm(niTz/L) sm(nirz'/L) f kpk p \ 

~ J (27r) d £ 2-/ ~ k 2 _ m 2 n + i£ ^ J (42j 

A corresponding result is obtained with Neumann boundary conditions. As before, we 
then drop the mode with n = 0. 



3.2 Electromagnetic expectation values 

We are now in the position of calculating the value of the energy-momentum tensor 
between the two plates. First we need the expectation value 



(Fp P (x)F^(x')) 



D 



nm'npu - Vauvpp)9 G D (x,x') 



D 



J \Vap,Vf3i> ~ VaPV0p)d z G D (x,x') 



(43) 



since the massless correlator satisfies the free wave equation (d 2 - d 2 )G D = 0. The 
structure of this result follows directly from antisymmetry of the field tensor and Lorentz 
invariance in the barred directions. As expected, it is simply given by the scalar correlator. 



In the same way we also find 

\Fp z (x)F Dz {x')) D = rjpp(d 2 - ^d 2 ^jG D (x,x') = — - \ ]p D d 2 z G N {x,x') 



D 



(44) 



The point-split expectation value of the Lagrangian density C = -F? v /A - F? z /2 is there- 
fore 



C(x,x')) D = -(D - 2)d 2 z G N (x,x') - G D (x,x') 



(45) 



in analogy with ([2 



From these results we can now read off the fluctations of the vacuum fields when the limit 
x — > x' is taken. For this purpose we combine (p7|) and (p8|) which give 



d 2 z G N/D = {D-l) 



( 47r )D/2 L D 



Cr(D) ± -f D (z/L) 



(46) 
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in this limit. With fi = v = in (I44p we then find for the z-component of the electric field 



\El) D = {D-2)- 



T{D/2) 



( 47r )D/2 L D 

The other components follow from f|43|) which gives 



(r(D) + -f D (z/L) 



, E i)D 



2T(D/2) 

\^)D/2 L D 



( R (D)--f D (z/L) 



(47) 



(48) 



where there is no summation over the index i^zon the left-hand side. For the fluctations 
in the magnetic components we similarly find (Bf z )o = — {E^)d and (Bfj)f) = — (Ef)£> 
where again there is no summation over the indices i, j ^ z. These relations also hold for 
the case of Neumann boundary conditions except for a change of signs in the last term of 
(JED and (gBD- 

These vacuum field fluctuations where first calculated by Liitken and Ravndal[6] for the 
ordinary Casimir effect with D = 4 spacetime dimensions in the Coloumb gauge and a 
different regularization. Using Cr(4) = vr 4 /90 and the function (I35p for f^z/L) in the 
above general results, we find 



{El) 



IT' 



48L 4 



(49) 



for Dirichlet boundary conditions, where 6 = ttz/L as before. The two other transverse 
components are given by (|48p as 



(El) 



48L 4 



F(0) - — 
K ' 15 



(50) 



Fluctuations of the transverse magnetic components are then (B£) = (By) = —{Ej) 

* z ) = -(El) 



while for the normal compoent we have (B^) = —(E^). 



Since the pressure is given by the expectation value of T z 
as 



F Z jiF^ z + C, it now follows 



P 



1 



(D - 2)d 2 z {G N + G D ) = (D-2)(D- 1)Eq 



(51) 



when expressed in terms of the energy density ([5]). It is again constant between the plates 
and a factor D — 2 times the scalar pressure (|26|) . This is exactly as expected since the 
Maxwell field has D — 2 scalar degrees of freedom. 

For the other components of the energy-momentum tensor we similarly find 



(Tp,p)D/N 



1D-2 
2D — 1 



-{D-2) 



G N + G D ±{D-A){G N -G D ) 
T{D/2) 



( 4vr )Z)/2 L D 



Cr(D) ± (D/2 - 2)f D (z/L) 



(52) 
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where the lower sign is for Neumann boundary conditions. It is only for D = 4 that the 
last, position-dependent term will be absent. And it is also then that Maxwell theory is 
conformally invariant. 

The above results for the electromagnetic field are very similar to what we found for the 
canonical, massless scalar field in the previous section. In that case the theory could 
be made conformally invariant with an improved energy-momentum tensor which gives a 
constant energy density. But for the Maxwell field there is no way to construct a local and 
gauge-invariant analogue of a similar Huggins term to cancel out the position-dependent 
part of (I52p . Thus the total Casimir energy obtained by integrating the energy-density 
is divergent and therefore looks different from what follows from the regularized sum of 
the zero-point energies of all modes. However, the difference turns out to be an infinite 
constant independent of the distance between plates. 

3.3 Discussion and conclusion 

The massless and free, canonical scalar field theory is not conformally invariant in other 
dimensions than D = 2. And it is only then that the Casimir energy density is constant and 
gives a finite integrated energy. The same satisfactory situation is also possible in higher 
dimensions when the theory is extended by making it conformally invariant, corresponding 
to adding the Huggins term to the energy-momentum tensor. 

This has been well-known for a long time, but not very well understood from a physical 
point of view. One of the most recent and detailed discussions of this phenomenon has been 
undertaken by Fulling who has attempted to understand the divergences in the canonical 
theory at a deeper level[17j. One can isolate the problem to the lack of commutativity 
between regularization of the integrated energy and the integration of the regulated energy 
density. This is perhaps not so surprising from a mathematical point of view, but hard 
to accept physically since the energy density is a physical quantity and should be tied up 
with the total energy of system. In other systems like the Casimir energy for a sphere, the 
energy density again diverges at the surface [18], but this is understood from its non-trivial 
geometry as first discussed by Deutsch and Candelas[19]. Since then the problem has been 
addressed by FullingfTT] and Milton with collaborators [20]. For a plane boundary there 
should be no such geometric complications. 

The electromagnetic Casimir effect for D = 4 is very similar to the scalar effect for D = 2. 
But for dimensions D > 4 there is no Huggins term for the electromagnetic case to cure the 
problem. From the point of view of the Casimir force alone, this is not a problem because 
the pressure is given by the expectation value ( T zz ) which is constant in all dimensions 
and equal to the force. But at first sight this force has little to do with the integral of the 
energy density (Too) which will always diverge at the plates when D > 4. 

This becomes especially clear when we just consider the electromagnetic fluctuations 
around one plate. The induced energy density can then be obtained from equation ([52]) 
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by taking the limit L — > oo. The pressure will then be zero on both sides of the plate 
while the other components become 

(T^d/n = TiP - 2)P/2 - 2) r( ^ } %p (53) 

since fr> = (L/z) D in this limit, as follows from the definition (|29p . It is non-zero on both 
sides of the plate and diverges when we approach it. The situation is analogous to the 
diverging energy density surrounding a classical pointlike electron. Thus, the behaviour 
([5*3")) is related to the intrinsic structure of a single plate, and the correponding integrated 
(infinite) energy is part of the energy required to make that plate. It does not contribute to 
the Casimir force. Thus, to find a connection between Casimir force and energy density it 
is sufficient to investigate the changes in energy density as two plates are brought together 
from infinite distance. We thus define 7^p as the expression (|52p subtracted contributions 
like (|53p from plates at z = (0, L), taking into account both sides of each plate. We find 
that 

rm/21 f T{D/2 ~ 2)(L/(L " z))D for z < °' 

r-^{z) = -(D - 2) \ * > w x { ( R (D)±(D/2-2)f D (z/L) for < z < L, (54) 
W L { T(D/2-2)(L/z) D for z > L, 

where fi)(z/L) = (h(D, 1 + z/L) + (h(D, 2 — z/L). This quantity T^p is finite everywhere, 
and its integrated energy agrees perfectly with the Casimir force since the integrals over 
the ^-dependent terms in (I54p cancel each other. The consistency between the various 
approches to the Casimir effect, i.e. total energy from mode sum ([2]), the pressure term 
(I51|) and the change in energy density (|54f) . gives support to the belief in a Casimir force 
which is essentially independent of the details of the plates. 

Of course, the interesting problem of the intrinsic and finite structure of a single plate re- 
mains. But this is similar to the problem of resolving the divergences caused by pointlike 
objects in quantum field theory. One obvious approach would be to modify the boundary 
conditions. If they were made softer so that they didn't affect fluctuations with wave- 
lengths below a certain cut-off A c , one would expect that the resulting energy density 
would be modified and finite for distances \z\ < X c away from the plate. This has actually 
been investigated by Graham et al. where a more physical boundary is described by an 
additional field [21]. 



In the more unphysical case where the Neumann or Dirichlet boundary conditions are 
replaced by periodic boundary conditions when D > 4, there would be no problems of 
these kinds. The energy density is then constant, giving a total energy consistent with 
the force between the plates. This is equivalent to the problem of photons in thermal 
equilibrium. Even if the trace of the energy-momentum tensor is non-zero when D > 4, 
the pressure P in this blackbody radiation is given by the energy density p by the standard 
expression P = p/d where d = D — 1 is the number of spatial dimensions [22J. 

But imposing such periodic boundary conditions, would be equivalent to just avoiding the 
problem. In conclusion, we must admit that the total vacuum fluctuations near confining 
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boundaries is still not completely understood, but is very likely to depend on the micro- 
scopic details of those boundaries. This is especially the case for the electromagnetic field 
in spacetimes with more than four dimensions. Fortunately the Casimir force seems to be 
rather insensitive to those details. 
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